The development of calibration and adaptive filtering procedures for neuromagnetic studies of human visual function

This thesis first considers the calibration and signal processing requirements of a neuromagnetometer for the measurement of human visual function. Gradiometer calibration using straight wire grids is examined and optimal grid configurations determined, given realistic constructional tolerances. Simulations show that for gradiometer balance of 1:104 and wire spacing error of 0.25mm the achievable calibration accuracy of gain is 0.3%, of position is 0.3mm and of orientation is 0.6°. Practical results with a 19-channel 2nd-order gradiometer based system exceed this performance. The real-time application of adaptive reference noise cancellation filtering to running-average evoked response data is examined. In the steady state, the filter can be assumed to be driven by a non-stationary step input arising at epoch boundaries. Based on empirical measures of this driving step an optimal progression for the filter time constant is proposed which improves upon fixed time constant filter performance. The incorporation of the time-derivatives of the reference channels was found to improve the performance of the adaptive filtering algorithm by 15-20% for unaveraged data, falling to 5% with averaging. The thesis concludes with a neuromagnetic investigation of evoked cortical responses to chromatic and luminance grating stimuli. The global magnetic field power of evoked responses to the onset of sinusoidal gratings was shown to have distinct chromatic and luminance sensitive components. Analysis of the results, using a single equivalent current dipole model, shows that these components arise from activity within two distinct cortical locations. Co-registration of the resulting current source localisations with MRI shows a chromatically responsive area lying along the midline within the calcarine fissure, possibly extending onto the lingual and cuneal gyri. It is postulated that this area is the human homologue of the primate cortical area V4

A peak-clustering method for MEG group analysis to minimise artefacts due to smoothness

Magnetoencephalography (MEG), a non-invasive technique for characterizing brain electrical activity, is gaining popularity as a tool for assessing group-level differences between experimental conditions. One method for assessing task-condition effects involves beamforming, where a weighted sum of field measurements is used to tune activity on a voxel-by-voxel basis. However, this method has been shown to produce inhomogeneous smoothness differences as a function of signal-to-noise across a volumetric image, which can then produce false positives at the group level. Here we describe a novel method for group-level analysis with MEG beamformer images that utilizes the peak locations within each participant's volumetric image to assess group-level effects. We compared our peak-clustering algorithm with SnPM using simulated data. We found that our method was immune to artefactual group effects that can arise as a result of inhomogeneous smoothness differences across a volumetric image. We also used our peak-clustering algorithm on experimental data and found that regions were identified that corresponded with task-related regions identified in the literature. These findings suggest that our technique is a robust method for group-level analysis with MEG beamformer images

A 3‐D Model of Gas Generation, Migration, and Gas Hydrate Formation at a Young Convergent Margin (Hikurangi Margin, New Zealand)

We present a three-dimensional gas hydrate systems model of the southern Hikurangi subduction margin in eastern New Zealand. The model integrates thermal and microbial gas generation, migration, and hydrate formation. Modeling these processes has improved the understanding of factors controlling hydrate distribution. Three spatial trends of concentrated hydrate occurrence are predicted. The first trend (I) is aligned with the principal deformation front in the overriding Australian plate. Concentrated hydrate deposits are predicted at or near the apexes of anticlines and to be mainly sourced from focused migration and recycling of microbial gas generated beneath the hydrate stability zone. A second predicted trend (II) is related to deformation in the subducting Pacific plate associated with former Mesozoic subduction beneath Gondwana and the modern Pacific-Australian plate boundary. This trend is enhanced by increased advection of thermogenic gas through permeable layers in the subducting plate and focused migration into the Neogene basin fill above Cretaceous-Paleogene structures. The third trend (III) follows the northern margin of the Hikurangi Channel and is related to the presence of buried strata of the Hikurangi Channel system. The predicted trends are consistent with pronounced seismic reflection anomalies related to free gas in the pore space and strength of the bottom-simulating reflection. However, only trend I is also associated with clear and widespread seismic indications of concentrated gas hydrate. Total predicted hydrate masses at the southern Hikurangi Margin are between 52,800 and 69,800 Mt. This equates to 3.4–4.5 Mt hydrate/km2, containing 6.33 × 108–8.38 × 108 m3/km2 of methane

Induced gamma activity in primary visual cortex is related to luminance and not color contrast: an MEG study

Gamma activity in the visual cortex has been reported in numerous EEG studies of coherent and illusory figures. A dominant theme of many such findings has been that temporal synchronization in the gamma band in response to these identifiable percepts is related to perceptual binding of the common features of the stimulus. In two recent studies using magnetoencephalography (MEG) and the beamformer analysis technique, we have shown that the magnitude of induced gamma activity in visual cortex is dependent upon independent stimulus features such as spatial frequency and contrast. In particular, we showed that induced gamma activity is maximal in response to gratings of 3 cycles per degree (3 cpd) of high luminance contrast. In this work, we set out to examine stimulus contrast further by using isoluminant red/green gratings that possess color but not luminance contrast using the same cohort of subjects. We found no induced gamma activity in V1 or visual cortex in response to the isoluminant gratings in these subjects who had previously shown strong induced gamma activity in V1 for luminance contrast gratings

Spherical harmonic based noise rejection and neuronal sampling with multi-axis OPMs

In this study we explore the interference rejection and spatial sampling properties of multi-axis Optically Pumped Magnetometer (OPM) data. We use both vector spherical harmonics and eigenspectra to quantify how well an array can separate neuronal signal from environmental interference while adequately sampling the entire cortex. We found that triaxial OPMs have superb noise rejection properties allowing for very high orders of interference (L=6) to be accounted for while minimally affecting the neural space (2dB attenuation for a 60-sensor triaxial system). We show that at least 11th order (143 spatial degrees of freedom) irregular solid harmonics or 95 eigenvectors of the lead field are needed to model the neural space for OPM data (regardless of number of axes measured). This can be adequately sampled with 75-100 equidistant triaxial sensors (225-300 channels) or 200 equidistant radial channels. In other words, ordering the same number of channels in triaxial (rather than purely radial) configuration may give significant advantages not only in terms of external noise rejection but also by minimizing cost, weight and cross-talk

SALT observations of the Chromospheric Activity of Transiting Planet Hosts: Mass Loss and Star Planet Interactions

We measured the chromospheric activity of the four hot Jupiter hosts WASP-43, WASP-51/HAT-P-30, WASP-72 & WASP-103 to search for anomalous values caused by the close-in companions. The Mount Wilson Ca II H&K S-index was calculated for each star using observations taken with the Robert Stobie Spectrograph at the Southern African Large Telescope. The activity level of WASP-43 is anomalously high relative to its age and falls among the highest values of all known main sequence stars. We found marginal evidence that the activity of WASP-103 is also higher than expected from the system age. We suggest that for WASP-43 and WASP-103 star-planet interactions (SPI) may enhance the Ca II H&K core emission. The activity levels of WASP-51/HAT-P-30 and WASP-72 are anomalously low, with the latter falling below the basal envelope for both main sequence and evolved stars. This can be attributed to circumstellar absorption due to planetary mass loss, though absorption in the ISM may contribute. A quarter of known short period planet hosts exhibit anomalously low activity levels, including systems with hot Jupiters and low mass companions. Since SPI can elevate and absorption can suppress the observed chromospheric activity of stars with close-in planets, their Ca II H&K activity levels are an unreliable age indicator. Systems where the activity is depressed by absorption from planetary mass loss are key targets for examining planet compositions through transmission spectroscopy

Retinotopic mapping of the primary visual cortex – a challenge for MEG imaging of the human cortex

Magnetoencephalography (MEG) can be used to reconstruct neuronal activity with high spatial and temporal resolution. However, this reconstruction problem is ill-posed, and requires the use of prior constraints in order to produce a unique solution. At present there are a multitude of inversion algorithms, each employing different assumptions, but one major problem when comparing the accuracy of these different approaches is that often the true underlying electrical state of the brain is unknown. In this study, we explore one paradigm, retinotopic mapping in the primary visual cortex (V1), for which the ground truth is known to a reasonable degree of accuracy, enabling the comparison of MEG source reconstructions with the true electrical state of the brain. Specifically, we attempted to localize, using a beanforming method, the induced responses in the visual cortex generated by a high contrast, retinotopically varying stimulus. Although well described in primate studies, it has been an open question whether the induced gamma power in humans due to high contrast gratings derives from V1 rather than the prestriate cortex (V2). We show that the beanformer source estimate in the gamma and theta bands does vary in a manner consistent with the known retinotopy of V1. However, these peak locations, although retinotopically organized, did not accurately localize to the cortical surface. We considered possible causes for this discrepancy and suggest that improved MEG/magnetic resonance imaging co-registration and the use of more accurate source models that take into account the spatial extent and shape of the active cortex may, in future, improve the accuracy of the source reconstructions

Can you tell your clunis from your cubitus? A benchmark for functional imaging

Advances in functional brain imaging have allowed the development of new investigative techniques with clinical application—ranging from presurgical mapping of eloquent cortex to identifying cortical regions involved in religious experiences. Similarly a variety of methods are available to referring physicians, ranging from metabolic measures such as functional magnetic resonance imaging and positron emission tomography to measurements based on electrical activity such as electroencephalography and magnetoencephalography. However, there are no universal benchmarks by which to judge between these methods. In this study we attempt to develop a standard for functional localisation, based on the known functional organisation of somatosensory cortex. Studies have shown spatially distinct sites of brain activity in response to stimulation of various body parts. Generally these studies have focused on areas with large cortical representations, such as the index finger and face. We tested the limits of magnetoencephalography source localisation by stimulation of body parts, namely the clunis and the cubitus, that map to proximal and relatively poorly represented regions of somatosensory cortex

The variance on the shock response of a carbon fibre composite due to the orientation of the weave

Three different orientations of a tape-wrapped carbon fibre composite with phenolic resin matrix (abbreviated to TWCP) have been investigated under one-dimensional shock loading. This has been achieved via a single-stage gas gun, with manganin gauges as the diagnostic tool. The orientations of TWCP studied in this paper were 25°, 45° and 90°, with respect to the impact face. The shock response of these orientations, for this material, has been obtained (the Hugoniot equation of state). These results have been contrasted with previously reported literature data for the same material at different orientations (0° and 20°). It was found that orientation had minimal effect on the behaviour of this composite under shock. The exception to this was the 90° orientation which exhibited an elastic precursor at particle velocities of less than 0.65 mm µs−1; where the shock velocity was equivalent to the elastic sound speed of the material

Ischaemia-reperfusion injury impairs tissue plasminogen activator release in man

AIMS: Ischaemia-reperfusion (IR) injury causes endothelium-dependent vasomotor dysfunction that can be prevented by ischaemic preconditioning. The effects of IR injury and preconditioning on endothelium-dependent tissue plasminogen activator (t-PA) release, an important mediator of endogenous fibrinolysis, remain unknown. METHODS AND RESULTS: Ischaemia-reperfusion injury (limb occlusion at 200 mmHg for 20 min) was induced in 22 healthy subjects. In 12 subjects, IR injury was preceded by local or remote ischaemic preconditioning (three 5 min episodes of ipsilateral or contralateral limb occlusion, respectively) or sham in a randomized, cross-over trial. Forearm blood flow (FBF) and endothelial t-PA release were assessed using venous occlusion plethysmography and venous blood sampling during intra-arterial infusion of acetylcholine (5-20 µg/min) or substance P (2-8 pmol/min). Acetylcholine and substance P caused dose-dependent increases in FBF (P<0.05 for all). Substance P caused a dose-dependent increase in t-PA release (P<0.05 for all). Acetylcholine and substanceP-mediated vasodilatation and substanceP-mediated t-PA release were impaired following IR injury (P<0.05 for all). Neither local nor remote ischaemic preconditioning protected against the impairment of substance P-mediated vasodilatation or t-PA release. CONCLUSION: Ischaemia-reperfusion injury induced substanceP-mediated, endothelium-dependent vasomotor and fibrinolytic dysfunction in man that could not be prevented by ischaemic preconditioning. CLINICAL TRIAL REGISTRATION INFORMATION: Reference number: NCT00789243, URL: http://clinicaltrials.gov/ct2/show/NCT00789243?term=NCT00789243andrank=1